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We report a microscopic and general theoretical 
formalism for electrical response which is appropri- 
ate for both DC and AC weakly nonlinear quantum 
transport. The formalism emphasizes the electron- 
electron interaction and maintains current conserva- 
tion and gauge invariance. It makes a formal con- 
nection between linear response and scattering ma- 
trix theory at the w^eakly nonlinear level. We derive 
the dynamic conductance and predict the nonlinear- 
nonequilibrium charge distribution. The definition of 
a nonlinear capacitance leads to a remarkable scaling 
relation w^hich can be measured to give microscopic 
information about a conductor. 
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Over the last decade tremendous effort have been de- 
voted in developing quantum transport theories which 
are appropriate for quantum coherent conductorstJ. 
While substantial progress have been achieved, so far 
there is still a lack of a general formalism which works 
not only in the linear DC regime but also in both DC and 
AC nonlinear regimes. Clearly, due to the great impor- 
tance to technological applications of quantum conduc- 
tors, such a formalism is urgently in need. AC transport 
involves time dependent fields, thus induction is impor- 
tant as characterized by, e.g. the existence of displace- 
ment current. Many transport theories do not consider 
this ingredient, resulting to a violation of current con- 
servation as observed by BiittikerEl. The nonlinear DC 
transport coefficients appear in front of powers of ex- 
ternal bias voltage in the expression of electric current: 
/ = I{{Va}) where Va is the potential at a probe labeled 
a which connects to the electron reservoir with chemi- 
cal potential /Iq. Thus a correct theory must maintain 
gauge invariance in addition to the current conservation: 
the physics should not change when potential everywhere 
is shifted by the same constant amount. This has severe 
constraints on the nonlinear DC transport coefficientscl. 

So far current-conserving and gauge invariant quan- 
tum transport has been analyzed using the scattering 
matrix Iheory (SMT) as developed by Biittiker and co- 
workergjD. Within this formalism, by including electron- 
electron (e-e) interaction, current conservation is ob- 
tained up to linear order of AC frequency lo for the dy- 
namic conductance Gapioj), and for special cases to the 
LiP term. The same interaction has led to a gauge in- 
variant second order weakly nonlinear DC conductance. 



In SMTcl, various ^-independent partial density of states 
(PDOS) appear naturally which characterizes the scat- 
tering. The SMT is quite intuitive and cas be imple- 
mented numerically for practical calculationsQ. However, 
so far no theoretical formalism exists which goes beyond 
the linear-w AC and second order DC transport coeffi- 
cients, and which satisfies current conservation and gauge 
invariance. There is also no approaches available to an- 
alyze the weakly nonlinear AC transport. Hence up to 
now these important problems have not been investigated 
systematically. 

It is well known that the linear coherent quantum 
transport can be discussedpusing either scattering ma- 
trix or using linear responseOij. In the linear DC rcgims- 
connection of these two approaches has been formalizedcl 
and attempt has been made to generalize it to nonlinear 
DC situationEl. Since a response theory is usually easier 
to be generalizable to higher order as it is in a perturba- 
tive formEl, we thus expect a general theoretical formalism 
derivable from the response theory which is appropriate 
for both DC and AC weakly nonlinear regimes in addi- 
tion to linear situations. We have indeed achieved this 
goal and it is the purpose of this Letter to report this 
development. 

Before go into the details, we summarize the main re- 
sults: (1) A general formalism for deriving AC and DC 
transport coefficients has been found which is current 
conserving and gauge invariant. (2) The formalism is 
put into a form which is numerically calculable for meso- 
scopic or even microscopic conductors order by order in 
weak bias and to all orders in lo. (3) The results obtained 
from SMT are contained in this formalism thus we make 
the formal connection between SMT and the response 
theory at the weakly nonlinear DC and AC level. (4) 
Generalized notions of the characteristic potential, the 
nonequilibrium charge, and the Lindhard functions are 
derived naturally from the time dependent internal re- 
sponse. (5) As examples of this formalism, we derived 
the linear (in bias) and second order weakly nonlinear 
dynamic conductance to all orders of w. We also predict 
the behavior of a nonlinear "capacitance" , and derive 
a remarkable scaling relation between the capacitances 
which is experimentally measurable and which gives mi- 
croscopic information such as the DOS of a system. 

A respon^ theory considers a time dependent 
perturbationD H' = V{r,t) = V{r)e-^\*\cos{nt) to the 
Hamiltonian Ho of a conductor. Here Ho contains the 
kinetic energy, the influence of an external static mag- 
netic field, the confining potential of the conductor, or 



any other static potential. We assume that the problem 
of Ho has been solved, Hgipn = ^n'4'm and we are in- 
terested in the effect of H'. It is crucial to realize that 
just solving the problem of Ho + H' will not generate a 
current conserving and gauge invariant transport theory, 
since the time varying field induces an internal potential 
U{{V}, t) which is a functional of the external perturba- 
tion V. It is this internal potential which generates such 
effect as a displacement current. Hence one must solve 
a quantum mechanical problem in conjunction with the 
Maxwell equations. We thus start from a self-consistent 
(up to the Hartree level) Hamiltonian H coupled with 
the Hclmhotz equation: 
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where Sp is the total charge involved in the response in- 
cluding those of the induced. 

We shall be interested in a multi-probe conductor small 
enough such that quantum coherence is maintained. In- 
side a probe a far away from the scattering region, the 
amplitude of the external disturbance is written as Vq. 
We shall calculate the internal potential C/ in a series 
form, expanded in powers of Va using the characteristic 
potentials (CP) Uoi{r,t), Uaf3(j,t) etc: 
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This expansion makes sense at weakly nonlinear regime. 
The CP characterizes the internal potential response to 
an external time-dependent perturbation. The physics 
associated with weakly nonlinear quantum transport is 
related to the CP's of the appropriate order (see below). 
Gauge invariance puts constraints on the CP, for instance 
X;„Ua(r,t) = 1, Y:^^Ua/3{r,t) = X;^Ma/3(r,i) = 0, and 
in general J2-y£8'^a{i3}i — 0. Here the subscript {P}i is a 
short notation of I indices 7, S,ri,- ■ ■. The quantum me- 
chanics problem of Eq. (y) is solved in the standard series 
fashion by iterating thci-Liouville-von Neumann equation 
for the density operators p(r, i) — pQ+J2i=iSpi{r,t). Us- 
ing the Liouvillc-von Neumann equation, order by order 
we thus derive equations for the terms Spi (r, t) . These 
equations can be formally solved in frequency space in 
terms of the characteristic potentials at the appropriate 
order: Spi is related to Uq,, Sp2 is related to Uap and Mq, 
5p3 to Wq/37 and the lower order w's, etc. Then using 
the density operator we can compute physical quantities 
such as..the charge density distribution and the electric 
currenija. 

With the charge density distribution and the Helmhotz 
equation (g), we can derive equations for the character- 
istic potential which is one of the central results of this 
work: 
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The right hand side of Eq. (0) is given by the frequency 
dependent local density of states (LDOS) which can be 
derived using the Green's functions. For example the 
lowest order (one index) LDOS is given by 
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whereQ W„ 



= '0n(r)E'V^m(i") with V the standard 
double-sided derivative operator; e„m = ^n — ^m, fnm = 
fi^n) — /(cm), and 77 is infinitesimal. The integration in 
(^) is along the boundary of probe a, and unit vector xi^ 
is the direction along this probe. This expression can be 
further written in terms of the Green's function. The 
third term on the left of Eq. (Q) is the induced charge 
written in terms of a frequency dependent Lindhard func- 
tion which is given by 
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Notice that the same Lindhard function appears in all 
equations. The boundary conditions can be set such that 
deep inside a reservoir connected to probe a, where the 
e-e interaction is completely screened, Ua — 1 while all 
others equal to zero. 

Some discussions of Eq. (Q) are needed. (1) The char- 
acteristic potential has multiple indices which are nec- 
essary in order to study nonlinear (in bias) transport 
coefficients. (2) All the terms are w-dependent as re- 
quired for a general theory which is appropriate to all 
orders in w. It is the combination of multiple indices 
and cj-dependence of the CP which allows the analy- 
sis of weakly nonlinear DC and AC transport in gen- 
eral terms of to. (3) The Lindhard function as derived 
is calculable from single particle states. It recovers the 
static result of Ref. || by letting lo ^ 0. For a perfect 
one-dimensional ballistic wire it is consistent with that 
of Ref. |l^ (including the double time derivative of Eq. 
(0)). (4) The concept of LDOS is generalized to higher 
order: at linear order (one index) it is the frequency de- 
pendent injectivityB, while at higher order (more than 
one indices) it also contains some internal response of 
the lower order. While their expressions (not shown ex- 
cept Eq. (^)) are complicated, they are all expressed in 



terms of the Green's functions. (5) It can be shown that 
the Lindhard function is related to the first order (one in- 
dex) total LDOS, J drill{r,ri,uj) — dn{r,uj)/de, where 
dn{r, u))/de = J2a driair, uj)/de. 

After solving the characteristic potentials of order / 
from Eq. (Q) , we thus obtain the total charge distribution 
at this order, 5pi[r,uj), which is just the right hand side 
of (|4|) subtracting the term with the Lindhard function, 
multiplying the external potential variations VaVfiV^ ■ ■■ 
(see Eq. (|^)). Various transport properties can be ob- 
tained immediately. 
Nonlinear "Capacitance". From the total charge 



distribution, we can 



l/2S/37^"/37^/5^7 



write Qa = Y^p'^apVp + 
■, where Qa is just the appro- 
priate spatial integration of the charge density. Hence 
the nonlinear theory naturally allows the definition of 
nonlinear capacitance coefficients. 
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dv 
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where index I > 1. The spatial integral over T means in- 
tegrating over the region where the charge Q is positive 
(or negative). At linear order in both to and voltage, i.e. 
when I = 2 andJetting w = 0, this gives the electro chemi- 
cal capacitancea which is of great experimental interesttil. 
Note that due to the finite screening length resulting from 
small DOS for mesoscopic conductors, a nonequilibrium 
charge distribution can be established even when there is 
DC coupling between the two capacitor "plates" (the +Q 
and -Q regions )t3. If we keep the general uj dependence 
but still work on the linear bias order {I — 2), Eq. (|^) 
gives the general linear dynamic response of the charge. 
A further very interesting result of ([7|) is the nonlinear 
"capacitance" for I > 2 and setting lu — Q. Such a quan- 
tity does not seem to have a geometrical counterpart, 
and it measures the degree of the nonequilibrium charge 
pile-up at the nonlinear order. 

In particular, let's calculate CafS-y for a parallel plate 
capacitor, where each plate has an area A and is in- 
finitely thin, and they are located in space on the y — z 
plane at positions x = and x = a. Using Eq. (0) 
to solve the characteristic potentials at different regions 
from X — — oo to x = +oo, matching the solutions at 
X = 0, a, we obtain ui and un. Within the Thomas- 
Fermi approximation of the Lindhard function, from Eq. 
(0) we obtain 
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where dNi/de is the total DOS on plate i, and Cn = 
[4:Tra/A+J2i=i \i{^'^d,Ni/de)]^^ is the usual electrochem- 
ical capacitancea. Eq. (H) is a remarkable relation be- 
cause on the left hand side there are only macroscopic 



quantities while on the right side all are microscopic. We 
have checked that for another system, being two large 
metallic rods with their ends at a distance a apart, ex- 
actly the same scaling of Cm/Cii is obtained, the only 
difference being that the factor 2 on the left side changes 
to 6. Hence this scaling relation gave us an extra handle 
on the microscopies of a conductor: by measuring capac- 
itances and forming the scaling combination, all that is 
left is the energy derivatives of DOS. For a symmetrical 
system Eq. (g) gives Cm — 0. This is because Cm is 
the coefficient of the term quadratic in bias, thus it must 
vanish for symmetric systems because the value of the 
pile- up charge cannot change when Va — > —Va for such 
a system. 

Linear Dynamic Conductance. As a second example 
we derive the dynamic conductance Gapi^uj). Gap is de- 
fined by the electric current flowing through the probe a: 
la = ^o Gafi{i^)Vfj. The current is obtained by a spa- 
tial integration of the current density across the trans- 
verse direction of the probe, with the current density 
given by Ji(r,w) = I]nm('^Pi('^))"'"(Jop(''))»™- Here 
(Jop(r))nm = —i{eTi/2m)'W nm- With these deflnitions, 
it is tedious but straightforward to evaluate Gajsiuj). The 
flnal result can be cast into a very compact form exactly: 



Gap{u:) = Gap{u: = 0) - iu:e^ / de(-/'(e)) 
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Here Gapii-o — 0) is the familiar linear DC conduc- 
tance whi«h can be calculated using the transmission 
coefficient^. g{Y^Yi,Lo) is Green's function of Eq. (|j), 
dna{vi,u})/de is another LDOS dual of dna{ri,u})lde^ 
and /'(e) — df /de. Setting to — Q, this reauJt reduces 
exactly to the emittance obtained by SMTB'EI. The lj- 
dependent parts are given as a sum of two contributions. 
The first is due to the external perturbation at a reser- 
voir, and it is determined by the w-dependent total PDOS 
dNaf3{oj)/de. The second is due to induction and is de- 
termined by the internal response. Although our theoret- 
ical formalism can be proven to conserve current, this is 
most easily seen from Eq. (O): it is not difficult to verify 
J^a^apii^) = 0. Finally it can also be shown that the 
following Onsager relation holds for the dynamic conduc- 
tance: Gapiui, B) = Gi3a{ui, — B). 

Weakly Nonlinear Dynamic Conductance. 

Our theory can be used to analyze weakly nonlinear AC 
transport to higher order in bias, and in the following we 
derive the second order expression defined by the second 



order electric current la 



(2) 



E, 



li^^ is 



calculated in similar fashion as the la of the last exam- 
ple, but using the density matrix in second order of bias 
The final result is 



GaPti^) — Gal3f{0J 



0)-i — UJ 
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de 
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de 
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In this result, GaPfit-^ = 0) is the second order weakly 
nonliaear DC conductance which was studied using 
SMTEI. Hence setting w = (|lO|) reduces to the known 
SMT result. ^^ '^-^'aB-y/'^^ equals the spatial integra- 
tion of the second order nonlinear LDOSEj dnp-y{v, uj)/d€. 
Result Eq. (nOh is general in frequency lo. The cur- 
rent conservation and gauge invariance can be explic- 
itly confirmed: Y.-y^aMii'^) = 0; I]aG"/57(^) = 
'^aGap-yioj) = 0. In addition, we have checked that 
there seems no simple relationship between GQ/3-y(w,B) 
and Gq/3^(w, — B), la agreement with the previous con- 
clusion on this pointu. 

Keeping terms linear in w, from Eq. ( |lO[) we de- 
rive the formula for the problem of second order Har- 
monic generation examined in Ref. |7| (which was done 
at w = 0). For the example of double-barrier resonant 
tunneling device, near a resonance energy E '^ Er, the 
scattering matrix takes the Breit-Wigner form. This 
allows simple expressions for the CPs: Ua — r^/r, 
Mil = -2e2(rir2/r2)(Je/|A|2), where Se = (E - Er), 
A = 6e + i(r/2), Fq, is the decay width due to barrier a, 
and F = Fi + F2. Evaluating all the terms of Eq. ( pO| ) 
using these expressions, we derive 



Giii(^) .^ F 

Giii(. = 0) - ^'^^2 



|A|2 



1 

f2 
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where Giii(w =|Q) is the second order weakly nonlinear 
DC conductancea. This result is also very interesting: 
it suggests that by measuring the ratio of the nonlinear 
conductances one obtains the microscopies of a tunneling 
device. 

The above three examples demonstrate the power of 
the present theoretical formalism: it is suitable for ana- 
lyzing weakly nonlinear DC and AC transport of meso- 
scopic conductors systematically. The expressions (M) 
and (llOl ) are general to all orders |0f uj. Setting w = 
they recover exactly those of SMTcl, thus providing the 
connection between SMT and the response theory at the 
weakly nonlinear AC level. In our formalism, the mul- 
tiple indexed characteristic potential describes the vari- 
ation of the scattering potential landscape of the con- 
ductor, at the appropriate weakly nonlinear level, as 
an external time dependent perturbation is applied to 
the electron reservoirs. The multiple indexed LDOS de- 
scribes the induced charge at the weakly nonlinear level, 
and they are related to the scattering Green's functions. 
For a megascopic conductorEI or even an atomic scale 
conductoiilil, the scattering Green's function can be eval- 
uated numerically thus our theoretical formalism pro- 
vides the basis of numerical analysis for a variety of quan- 
tum conductors. Finally, we comment that in this pa- 
per we have used the Lorentz gauge for electrodynamics, 
we have made an explicit check that using the Coulomb 
gauge gives exactly the same result. Clearly any gauge 
should work as is required by the Maxwell equations. 



Central to the theory presented here is the self- 
consistent coupling of the quantum mechanical equation 
with the Maxwell equations, thus the internal response of 
the system is taking into account which is crucial to ob- 
tain electric current conservation and gauge invariance. 
Conceptually, we have introduced the multiple indexed, 
frequency dependent characteristic potential and the lo- 
cal density of states. We have also extended the concept 
of nonequilibrium charge distribution to the nonlinear or- 
der. Our response theory is appropriate to equilibrium 
or near-equilibrium properties. For far-from-equilibrium 
problemS|-one may employ the nonequilibrium Green's 
functionstJ. A simple application of our theory naturally 
leads to the concept of nonlinear capacitance. A remark- 
able scaling relation is predicted such that a macroscopi- 
cally measurable ratio of these capacitances is related to 
the microscopic information of a conductor. This result 
suggests an interesting experiment by using two small 
conductors which couple capacitively, and measure the 
AC response as a function of the amplitude of the AC 
bias. 
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